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昆虫 种 群 累积 发 育 频数 的 数学 
模拟 中 时 沸 个 体 的 模拟 处 理 


《北京 农业 大 学 》《 中 国 科学 院 动 物 研 究 所 》 


昆虫 种 群生 态 过 程 的 数学 模拟 一般 多 以 研究 单 种 种 群 在 某 一 世代 不 同 发 育 阶 段 中 
的 数量 变动 为 基础 。 一 个 种 群 在 某 一 世代 中 从 某 一 发 育 阶 段 进 人 到 另 一 发 育 阶段 时 , 由 
于 个 体 发 育 速率 的 差异 ， 形 成 了 一 种 发 育 时 间 的 分 布 。 de Wit 等 (1974)、Goudriaan 
《1973)、Piclou 《1969) 和 Stinner 等 (1975) 曾 先后 提出 几 种 拟 合 发 育 时 间 分 布 曲线 的 方 
法 ， 但 都 没有 对 种 辩 内 由 于 小 部 分 上 发育 特别 退 缓 的 个 体 所 造成 的 发 育 时 间 分 布 的 显著 偏 
态 问题 进行 讨论 。 

本 文 以 棉铃 虫 在 不 同 温度 下 各 发 育 期 的 实验 室 资 料 为 基础 ， 提 出 一 稳 数 字模 拟 过 程 
中 处 理发 育 迟 缓 个 体 (简称 时 滞 个 体 ) 的 方法 。 

一 、 间 题 的 提出 

一 个 种 群 的 增长 过 程 ,是 种 群 中 每 个 个 体 生长 发 育 状 态 的 总 表现 ,由 于 个 体 之 间 的 差 
异 , 发 育 有 先 有 后 ,发 育 较 早 和 较 晚 的 个 体 少 ,中 间 类 型 的 较 多 , 形成 “两 头 少 , 中 间 多 ”的 
普遍 规律 ,根据 视 铃 虫 在 不 同 温度 不 同 发 育 阶段 的 资料 分 析 , 还 表明 种 群 中 的 “两 头 ” 并 不 
对 称 , 大 多 数 表 现 为 发 育 时 间 偏 长 的 个 体 数 多 于 发 育 时 间 偏 短 的 个 体 数 。 往 往 由 于 存在 一 
小 部 分 个 体 发 育 特别 迟缓 ,把 整个 种 群 的 发 育 时 间 分 布 拉 得 很 长 的 现象 ( 见 表 1)。 对 于 这 





表 1 在 邱 "C 条 件 下 ,棉铃 虫 幼虫 化 晴 时 间 的 分 布 ( 引 自 中 国 科 学 院 动 物 研究 所 昆虫 生态 室 实 验 ) 


化 旺 时 间 ( 天 ) 频 数 《个 ) 


16 一 17 
17 一 18 
18 一 19 
19. 一 20 
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种 发 育 时 间 的 分 布 资料 ,我 们 曾 采 用 Goudriaan 《1973) 和 de Wit 等 (1974) 的 对 称 分 布 
来 拟 合 ( 见 表 2) ， 模 拟 的 结果 经 过 xX- 检验 ,不 应 认为 表 1 资料 具有 Goudriaan《1973) 和 
de Wi 等 (1974) 所 给 出 的 对 称 分 布 的 形式 。 我 们 亦 用 Sinner《〈19757 提出 的 偏 态 分 布 对 
表 1 资料 进行 拟 合 , 仍 出 现 “ 两 头 ”理论 个 数 偏 多 ,“ 中 间 ” 理 论 个 数 偏 少 的 现象 ,结果 也 不 
理想 ( 见 表 3)。 因 此 , 有 必要 提出 一 种 适用 于 模拟 如 表 1 所 给 出 的 资料 的 分 布 模型 。 

训 2 根据 Goudriann (1973) 的 方法 对 表 1 资料 的 模拟 结果 的 统计 分 析 


区 
O 沪 
酒 


天 数 
《组 中 值 》 







15.5 -0 } 6.0036 
16.5 9.2 
17.5 18.5 1.6351 
18.5 21.3 2.7836 
19.5 16.7 0.6521 
20.5 9.6 4.5375 
21.5 4.4 
22.5 1.6 
23.5 0.5 0.4313 
24.5 0.2 
25.5 0.0 


自由 度 为 6 一 2 一 1 一 3 的 和 变量 在 信 度 为 0.05 时 的 检验 临界 值 总 一 7.815。 故 小 于 3 (2 二 祁 ， 应 认为 
理论 分 布 与 实际 分 布 差异 显著 


朗 3 根据 Stinner 等 (1975) 提出 的 模型 对 训 1 资料 的 氛 合 结果 的 统计 分 析 
1 理论 累积 频数 
84.y， 








12. 4 12.44 3 7.1635 

29.30 16.86 24 3.0237 

45.99 16.69 29 9.0795 

060.02 14.03 20 2.5403 

70.43 10.41 5.2746 
5. 





er 
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:33.1919 





自由 度 为 7 一 2 一 1 一 4 的 入 变量 在 信和 度 为 0.05 时 的 检验 临界 值 X5 一 9.488。 故 小 于 互 2 应 认为 


理论 分 布 与 实际 分 布 差 异 显著 。 
中 =， 一 3 二 二 二 (1 一 5s) 其 中 B 一 235.5，4 一 15.5，9 =1.02523, 及 二 2.01137。 





二 、 假 设 与 模型 
表 1 的 资料 表明 ,该 种 群 中 大 多 数 个 体 在 第 21 天 以 前 化 师 , 在 第 18 一 19 天 之 间 化 晴 
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的 个 体 数 最 多 ;而 且 大 多 数 个 体 的 发 育 时 间 对 称 地 分 布 在 这 个 众 数 时 间 的 附近 ( 兄 表 工 虚 
线 上 部 分 )。 另 外 ,有 一 小 部 分 个 体 发 育 特别 迟缓 ; 零 零星 星 地 拖 到 第 26 天 才 完 全 化 晴 。 
为 讨论 方便 ,把 这 些 个 体 叫 做 时 灌 个 体 ”。 

显然 .一 个 种 群 的 动态 主要 决定 于 大 多 数 个 体 , 表 1 中 占 总 数 只 6% 左右 的 时 灌 个 体 
对 整个 种 群 的 动态 是 不 起 决定 性 作用 的 。 ee 
分 开 , 将 整个 种 群 看 成 是 两 个 总 体 的 组 合 : 一 个 是 由 正常 个 体 组 成 的 总 体 ,它们 是 影响 种 
群 动态 的 决定 性 因素 ; 另 一 个 是 由 时 海 个 体 所 组 成 的 总 体 ， 它们 只 能 对 整个 种 群 的 动态 起 
修正 作用 。 例 如 表 1 中 大 多 数 个 体 在 最 多 一 天 ( 即 第 18.5 天 ) 的 前 后 各 两 天 时 间 内 化 师 ， 
也 可 以 说 在 第 21 天 以 前 化 晴 , 故 可 把 它们 当成 一 个 正常 总 体 ;而 化 里 时 间 在 21 天 以 后 的 
个 体 ,就 划 归 时 兹 总 体 。 这 样 划分 后 ,基本 上 就 可 以 把 正常 总 体 当 成 具有 对 称 分 布 形式 的 
了 。 

为 检验 这 种 划分 的 合理 性 ,可 将 正常 总 体 的 发 育 时 间 期 望 值 及 均 方 差 的 估计 值 求 出 ， 


分 别 得 到 六 二 二 18.45， = I $ 二 0.9323, 故 正常 总 体 的 发 育 时 间 应 该 在 次 十 
了 一 


36 一 21.25 (天 ) 以 内 。 可 见 化 旺 时 间 在 第 21 天 以 后 的 个 体 , 不 应 与 正常 个 体 同 属 一 个 发 
育 时 间 分 布 总 体 之 内 **, 也 就 是 说 , 把 种 群 按 正 常 总 体 与 时 灌 总 体 划分 的 这 种 做 法 有 统 
计 检 验 作为 依据 。 

把 单 种 种 群 分 成 正常 总 体 与 时 潇 总 体 进行 处 理 , 就 是 本 文 提出 的 模型 的 基本 假定 。 下 
面 讨论 分 为 两 个 总 体 后 的 种 群发 育 时 间 分 布 模型 。 

假定 正常 总 体 与 时 滞 总 体 的 发 育 时 间 分 布 频率 函数 分 别 为 he) 及 f(z) 《其 中 (2) 
可 以 是 具有 对 称 分 布 形 式 的 分 布 频率 函数 ) 两 个 总 体 在 种 群 中 所 占 的 百分率 分 别 为 九 及 
和 《其 中 为 十 为 一 1)。 那么 ,正常 总 体 及 时 谐 总 体 的 发 育 时 间 分 布 浮 数 分 别 为 : 


FOOD 一 | farR PCD 一 | Mr)ar， G) 


若 以 F(z) 表示 整个 种 群 的 发 育 时 间 分 布 函 数 , 则 可 以 证 明 :; 
F(#) 一 和 PiGz) + hF,(z) 


一 | [ft) + hflr) lar, (2) 


从 而 得 到 种 群 的 发 谊 时 间 分 布 频率 项 数 为 : 
fa = Hf) 十 1j(9) 
= f(t) 十 22[ 户 oo — f(a)1. 
事实 上 , 当 种 群 的 总 个 体 数 六 很 大 , 且 正常 总 体 和 时 澡 总 体 的 个 体 数 Ni 及 N; 也 都 充 
分 大 时 ,可 以 根据 (1) 得 到 时 刻 * 种 群 中 正常 总 体 的 已 发 育 个 体 数 为 NF1(z)， 时 沸 总 体 
的 已 发 育 个 体 数 为 ;F(z)。 故 种 群 在 时 刻 * 的 已 发 育 个 体 数 为 NiF1C2) 十 N2F2《z), 这 
个 数目 应 该 等 于 种 群 的 总 数 久 乘 以 种 群 的 发 育 时 间 分 布 函 数 在 时 刻 t 的 值 F(z), 即 





* 种 群 中 存在 部 分 发 育 特别 迟缓 的 个 体 这 种 规律 狂 观 象 叫做 时 淆 现象 。 本 文 招 这 种 个 体 叫做 时 满 个 体 。 并 为 区 
别 起 见 , 把 非 时 兴 个 体 叫 做 正常 个 体 。 

* 在 21 一 22 天 中 化 师 的 个 体 ,应 认为 是 在 31.5 天 (>>21.25 天 ) 化 到 , 故 也 不 属于 正常 总 体 。 此 外 ,以 砍 十 36 
做 为 检验 的 界限 是 券 虑 到 正常 总 体 具 有 Goudriaan (1973》 所 给 出 的 对 称 分 布 这 一 假定 。 
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TIME Y Yanin Ymar 
0.00000 0.79002,。 十 2 

15.00 0.00000 :  w* 

15.50 0.00000 : 水 

16.00 0.00000 : 关 

16.50 0.28379,,+1 : 水 

17.00 0.11922,o 十 2 ; 六 

17.50 0.26788, 十 2 : * 

18.00 0.43379,o 十 2 : 米 

18.50 0.57580,o 十 2 : 水 

19.00 0.67526,o 十 2 : 来 

19.50 0.73462, 十 2 : 3 

20.00 0.76565,o 十 2 : * 

20.50 0.78016u 十 2 “ 水 

21.00 0.78632, 十 2 : * 

21.50 0.78872 十 2 : 米 

22.00 0.78959 +2 : 

22.50 0.78989,o 十 2 。 @ 水 

23.00 0.78998,。 十 2 。 来 

TIME Y Ya Ya 
0.00000 0.84000,o 十 2 

15.00 0.00000 : * 

15.50 0.00000 : 米 

16.00 0.00000 :水 

16.50 0.28379 十 1 。 * 

17.00 0.11922,,+2 ° 水 

17.50 0.26788,o 十 2 水 

18.00 0.43379o 十 2 : * 

18.50 0.57580,,+2 : 水 

19.00 0.67526,u 十 2 : 水 

19.50 0.73474,o 十 2 : 来 

20.00 0.76645,o 十 2 ， * 

20.50 0.78291o 十 2 : * 

21.00 0.79289o 十 2 : * 

21.50 0.79529,o 十 2 : * 

22.00 0.80188,。 十 2 : * 

22.50 0.80920.0+2 : * 

23.00 0.81665,o 十 2 : 水 

23.50 0.82342o 十 2 . 水 

24.00 0.82896,o 十 2 : 水 

24.50 0.83309 十 2 : * 

25.00 0.83309o 十 2 : * 

25.50 0.83591o 十 2 ， 水 

26.00 0.83771,o+2 ; ® 水 

26.50 0.83878,。 十 2 ， 六 

27.00 0.83938 +2 ， 六 


图 1 资料 的 模拟 结果 图 
(加 杀 指 C2) 式 所 给 模型 的 正常 总 体 的 模拟 输出 ,加 系 指 (2) 式 所 给 模型 的 整个 种 群 的 模拟 输出 》 


2 期 陈 维 博 等 : 昆虫 种 群 累积 发 育 频数 的 数学 模拟 中 时 灌 个 体 的 模拟 处 理 153 


NF(#) = NF() + NsF,(t), 
考虑 到 和 一 Ni/N 及 5 二 NWN , 则 有 : 
| F(2) = 人 PiCz) 十 hb,F(2). 
根据 (2) 不 难得 到 , 两 个 总 体 的 期 望 值 m; 及 方差 oii 一 1, 2) 与 整个 种 群 的 期 望 值 
妈 及 方差 中 之 间 存 在 下 列 关系 式 : 


m = hm 十 bm, (3) 
= ho? + ho + hm — m2), (4) 

(2) 式 可 进一步 化 为 : 
F(z) 一 F(#) + hl F(t) 一 F(2)]. (5) 


可 见 (2) 式 所 给 出 的 种 群发 育 时 间 分 布 模型 F(z), 实际 上 是 以 正常 总 体 的 对 称 分 布 已 (z) 
为 基础 ,并 附 有 一 个 修正 项 加 LF2tz) 一 F(z)]， 修 正 作用 的 大 小 由 为 值 的 大 小 决定 。 当 
久 值 比较 小 ,即时 滞 总 体 在 种 群 中 所 占 的 比例 很 小 时 ,种 群 的 发 育 时 闻 分 布 将 以 正常 总 体 
的 发 育 时 间 分 布 F(z) 为 主体 (参见 图 1 中 @ 与 @ 的 比较 )。 有 了 时 时 灌 个 体 只 有 一 、 两 个 ， 
此 时 可 以 忽略 修正 项 ,同时 种 群 分 布 也 以 对 称 分 布 的 形式 出 现 。 

表 4 根据 (2) 式 的 楼 型 对 表 1 资料 的 模 氛 结果 的 统计 分 析 

















07) (9) 6) 
理论 分 布 | 。 经 验 分 布 


4) 3 
F= (49 F* 一 一 
元 避 士 0) 一 六 









(2) 污 蒜 果 加 六 站 的 到 | G7 (© 
天 数 | 实际 需 生 PS% 论 累加 频 | 纹 放 开 数 [GDC] 
大 Tt T 名 pe NREC 一 1) (5) 






1FC2) 一 
FxCs + 0)| 


































15.5 0.0 0-0000 0.0000 0.0000 
16.5 2.8 0.01429 0.0333 0.0357 0.002+ 
17.5 26-8 0.00000 0.3190 0.3214 0.002+ 
18.5 57.6 0.10519 0.6857 0.6667 0.0190 
19.5 73.5 1.05723 0.8750 0.9048 0.0298 
20.5 78.3 0.67500 0.9321 0.9405 0.0084 
21.5 79.53 0.03333 0.9464 0.9524 0.0060 
22.5 80.9 0.11429 0.9631 0.9643 0.0012 
23.5 82.3 0.11429 0.9798 0.9762 0.0036 
24.5 83.3 1.00000 0.9917 0.9881 0.0036 
25.5 83.6 1.63333 0.9952 1.0000 0.0048 











~ 一 


3.74695 DN=max| F(2) ~ F*(14+0)| =0.0298 
了 


根据 对 棉铃 虫 在 不 同 温度 条 件 下 种 群 在 各 个 发 育 期 的 动态 资料 的 分 析 ， 可 以 假定 十 
总 体 的 绝 大 多 数 的 正常 总 体 , 符 合 Goudriaan 《1973》 提出 的 对 称 型 分 布 。 根 据 这 一 假定 
为 前 提 , 可 把 不 对 称 的 剩余 部 分 划 为 时 请 总 体 。 我 们 依照 上 述 假定 及 模型 ,在 处 理 时 汪 总 
体 的 分 布 函数 形式 时 , 近似 地 仍 采 用 与 正常 总 体 同样 的 对 称 型 分 布 。 从 不 同 温度 条 件 下 
各 个 发 育 期 的 模拟 结果 来 看 ,(2) 式 给 出 的 模型 与 实际 资料 相当 吻合 。 


三 、 结 有 果 和 讨论 z 
对 表 1 的 资料 以 21 天 为 界 划分 种 群 为 两 个 总 体 , 并 分 别 用 Goudriaan 1973) 所 提出 
的 模型 进行 模拟 ,再 按 (2) 的 模式 屋 加 两 个 总 体 ,作为 整个 种 群发 育 动态 的 模拟 ,得 到 了 比 
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较 满意 的 结果 。 1 是 正常 总 体 与 整个 种 群 的 电子 计算 机 模拟 输 出 结果 NF(z) 以 及 
NF(z)。 表 4 根据 表 1 的 实际 资料 , 与 图 1@ 的 模拟 数值 结果 做 比较 。 可 以 看 到 ， 理论 
分 布 F(z) 与 经 验 分 布 F*(z 十 0) 十 分 吻合 ( 见 表 4(3)、(4) 栏 )。 表 4 第 (6) 栏 进一步 根 
据 第 (2) 及 第 (5) 栏 求 出 各 时 间 区 间 对 应 的 冶 值 ,得 到 的 关 值 都 很 小 ,它们 的 总 和 也 相当 
小 。 由 于 (5) 栏 中 理论 频数 小 于 5 的 需 做 适当 的 合并 ,合并 后 只 有 四 项 ,但 决定 模型 C2) 的 
理论 分 布 的 参数 共有 入 , mi my, of 及 有 等 五 个 (参见 (3)、(4) 和 (5) 式 ), 参 数 个 数 多 于 
项 数 不 适 于 用 涛 值 做 分 布 适合 度 的 检验 。 这 里 进一步 采用 柯 尔 莫 哥 洛 夫 的 分 布 函数 拟 
合 优 度 检验 标准 来 检验 。 由 表 4 第 (4) 栏 得 到 : 

Dv 二 max| F(z2) 一 Fr 十 0)| = 0.0298 
也 柯 氏 检 验 表 可 查 得 , ”一 85，a 一 0.05 时 有 : 

Ds,,s = 0.13072 > 0.0298 = Dwo 
完全 可 以 认为 经 验 分 布 具有 理论 分 布 (2) 的 形式 。 

图 2 给 出 30%C 条 件 下 棉铃 虫 化 师 时 间 的 资料 及 采 衣 模型 (2) 的 模拟 结果 的 比较 ,可 

以 看 到 理论 值 与 实际 值 也 比较 吻合 。 


TIME Y YMIN YMAX 
0.00000 0.84000,o 十 2 

10.00 0.00000 : 来 

10.50 0.00000 : 图 

11.00 0.00000 : A* 

11.50 0.84061o 十 1 : 口水 

12.00 0.26936,, 十 2 : * 

12.50 0.4792410 十 2 : 口 水 

13.00 0.64241,, 十 2 : 水 

13.50 0.74079,o 十 2 米 口 

14.00 0.78905,0+2 米 

14.50 0.81099,0+2 米 口 

15.00 0.82106,。 十 2 * 

15.50 0.82669,o+2 图 

16.00 0.83066,o 十 2 水 

16.50 0.83373,0 十 2 胃 

17.00 0.83603,。 十 2 六 

17.50 0.83763,0 十 2 力 

14.00 0.83764,0 十 2 水 

18.540 0.83867,0 十 2 : 水 

19.00 0.83920,o 十 2 : 口 实际 水 

19.50 0.83964.0+2 : 炒 ” 模 所 * 

20.00 0.83983,o 十 2 : 困 ”实际 与 模拟 数据 相 结 合 米 


图 2 在 30%C 条 件 下 棉铃 虫 幼 虫 化 师 时 间 的 实际 分 布 及 模拟 结果 
根据 对 棉铃 虫 在 不 同 温度 下 不 同 发 育 期 的 模拟 分 析 , 可 以 认为 ; 
1. 存在 时 滞 现 象 的 昆虫 种 群 的 模拟 模型 ， 可 以 用 正常 个 体 构 成 的 总 体 与 时 沾 个 体 构 
成 的 总 体 组 合 而 成 的 群体 模型 。 在 固定 温度 条 件 下 , 种 群 的 发 育 时 间 分 布 可 以 用 两 个 总 
体 的 发 育 时 则 分 布 的 线性 组 合 表示 ,线性 组 合 的 系数 是 两 个 总 体 在 促 群 中 所 占 的 比例 ,在 
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做 适当 的 划分 后 ,两 个 总 体 的 发 育 时 间 分 布 都 可 采用 Goudriaan (1973) 所 提出 的 对 称 分 
布 形式 进行 模拟 。 

2. 在 变温 条 件 下 的 模拟 可 做 如 下 处 理 : 把 不 同 温 度 下 两 个 总 体 的 发 育 率 以 及 所 占 的 
比例 均 视 为 温度 的 函数 。 在 模拟 过 程 中 的 每 个 时 间 步 长 内 , 可 近似 地 把 发 育 率 和 两 总 体 
的 比例 当成 是 常数 《因为 这 时 温度 也 可 近似 地 当成 是 不 变 的 )， 这 样 就 可 如 de Wit 等 
《19747 的 方法 采用 线性 组 合 的 形式 来 处 理 。 

3. 鉴于 本 文 用 以 研究 分 析 的 资料 是 实验 室 的 资料 ,种 群 的 总 个 体 数 不 够 大 (一 般 只 在 
100 个 左右 )， 抽 得 结果 的 统计 分 析 受 到 一 定 的 影响 (例如 不 可 能 都 同样 用 六 检验 法 作为 
检验 分 布 拟 合 程度 的 统一 标准 ), 有 可 能 存在 某 些 偶然 误差 ,有 待 进一步 检验 。 虽 然 如 此 ， 
从 表 2、3、4 中 如 值 栏 中 可 以 看 到 ,各 相同 时 间 区 间 中 的 妈 值 除 个 别 外 :一般 表 4 中 的 值 
都 比 表 2 或 3 中 的 各 对 应 值 小 ,说 明 本 文 所 提出 的 模型 (2) 所 做 的 改进 确 是 可 取 的 。 
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A MATHEMATICAL SIMULATION TECHNIQUE FOR THE 
TIME-LAG INDIVIDUALS IN AN INSECT POPULATION 
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This paper deals with the digital simulation in different developmental stages 
of the eotton bollworn Heliothis armigera, at various temperatures. A simulation 
model, based on the study of the time-lag individuals in eaeh developmental stage, 
was proposed. 


